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ABSTRACT 

We present the discovery of radio afterglow emission from the high redshift (z = 6.295) burst 
GRB 050904. The peak flux density for this burst is similar to typical low-redshift gamma-ray bursts 
(GRB). We further show that beyond a redshift of order unity, the flux density of radio afterglows are 
largely insensitive to redshift, consistent with predictions. By combining the existing X-ray, near- infrared 
and radio measurements, we derive estimates for the kinetic energy and opening angle of the blast wave, 
and for the density of the circumburst medium into which it expands. Both the kinetic and radiated 
energy indicate that GRB 050904 was an unusally energetic burst (10 52 erg). More importantly, we are 
able to make an in situ measurement of the density structure of the circumburst medium. We conclude 
that GRB 050904 exploded into a constant density medium with n o =680 cm" 3 , which is two orders of 
magnitude above the nominal value for low-redshift GRBs. The next generation of centimeter (EVLA) 
and millimeter radio instuments (ALMA) will be able to routinely detect events like GRB 050904 and 
use them to study magnetic fields, and the atomic and molecular gas in the high redshift Universe. 

Subject headings: gamma-ray bursts: specific (GRB 050904) 



1. INTRODUCTION 

Understanding the reionization of the Universe, when 
the first luminous sources were formed, is one of the lastest 
frontiers of observational cosmology. Constraints have 
been obtained using diagnostics such as quasar studies of 
the Gunn-Peterson absorption trough, the luminosity evo- 
lution of Lya galaxies, and the polarization isotropy of 
the cosmic microwave background. Taken together, these 
data portray a complicated picture in which reionization 
has taken place over a wide range of redshifts (6< z < 20) 
rather than at one specific epoch. The dominant source 
of reionization appears to be due to ultraviolet emission 
from young, massive stars (see review by Fan, Carilli & 
Keating 2006). 

As the most luminous explosions in the Universe, 
gamma-ray bursts (GRBs) are potential signposts of these 
early massive stars. The radio, infrared and X-ray after- 
glow emission from GRBs are in principle observable out 
to z~ 30 (Miralda-Escude 1998; Lamb & Reichart 2000; 
Ciardi & Loeb 2000; Gou et al. 2004; Ioka & Meszaros 
2005). ft is estimated that 10% of GRBs detected by 
the Swift satellite are at z > 5 (Natarajan et al. 2005; 
Bromm & Loeb 2006). The most distant GRB to date is 
GRB 050904 at z = 6.295. It was detected by the Swift 



Burst Alert Telescope and localized to an accuracy of a 
few arcseconds by the Swift X-ray Telescope (Gehrels et 
al. 2004). Follow-up optical and near-infrared observa- 
tions (NIR) were begun shortly thereafter. Details on the 
discovery and properties of the X-ray and NIR afterglow 
of this burst can be found in several papers (Cusumano et 
al. 2006; Haislip et al. 2006; Boer et al. 2006). An opti- 
cal/NIR spectrum taken by the Subaru telescope 3.4 days 
after the burst showed multiple heavy metal absorption 
lines and a Gunn-Peterson trough with a damping wing 
redward of the Lya cutoff which was used to derive the 
neutral fraction xr < 0.6 (Totani et al. 2005; Kawai et 
al. 2006). With a Swift GRB detection rate of 100 yr" 1 , 
GRBs could one day replace quasars as the preferred probe 
of the high redshift Universe. 

In this paper we report on the detection of the radio af- 
terglow from GRB 050904 (§2), which makes it possible to 
derive physical properties of the explosion and the circum- 
burst medium (§3). Our results are compared with pre- 
dictions for the properties of the explosion and how the 
progenitors of high redshift GRBs are expected to have 
shaped their surrounding environments (§4). 

2. OBSERVATIONS AND RESULTS 
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Radio observations were undertaken with the Very 
Large Array 12 (VLA) at a frequency of 8.46 GHz (see 
Table 1). To maximize sensitivity, the full VLA contin- 
uum bandwidth (100 MHz) was recorded in two 50 MHz 
bands. Data reduction was carried out following standard 
practice in the AIPS software package. Some additional 
phase and amplitude self-calibration was necessary to re- 
move the contaminating effects of a bright (10 mJy) radio 
source four arcminutes northeast from the GRB. 

Radio emission from GRB 050904 was not detected dur- 
ing the first week, ruling out a bright short-lived compo- 
nent similar to GRB 990123 (Kulkarni et al. 1999). By av- 
eraging the data from this period, we obtain a peak flux at 
the position of the NIR afterglow of 33±14 uJy. A clear de- 
tection was made during three epochs 34-37 days after the 
burst (see Table 1). Averaging the data taken in October 
we obtain a mean flux density of 76±14 uJy. The spectral 
radio luminosity, expressed as L v — AirF^ d 2 L (1 + z) a_/3_1 
= 2.5xl0 31 erg s _1 Hz -1 (where F v ex £ C V 3 and a ~ and 
f3 = 1/3 has been assumed, corresponding to an optically 
thin, flat post-jet break light curve), is normal for GRBs 
(Soderberg et al. 2004), and is two orders of magnitude 
below the highest redshift radio-loud quasars (Carilli et al. 
2004). 

We show these detections in Fig. 1 together with a com- 
plete sample of 8.5 GHz flux density measurements for 
GRBs with known redshifts. Neither the time-to-peak nor 
the flux density of GRB 050904 is unusual compared to the 
known sample of GRBs at lower redshifts. In the source 
rest frame the flux density for the majority of afterglows 
is reached before 5 days postburst. More interestingly, the 
average centimeter flux density in Fig. 1 shows only a weak 
dependence on redshift. 

This effect was predicted by Ciardi & Loeb (2000) and 
is shown here for the first time. It is reminiscent of the 
"negative k-correction" for submillimctcr observations of 
ultraluminous infrared galaxies (Blain et al. 2002). The 
afterglow flux density remains high because of the dual ef- 
fects of spectral and temporal redshift, offsetting the dim- 
ming due to the increase in distance (Lamb & Rcichart 
2000). We illustrate in Fig. 1 how slowly the centimeter 
flux density decreases for a canonical GRB afterglow be- 
yond z ~ 1 due to this effect (long dashed lines) , compared 
to one whose luminosity is assumed constant (short dashed 
lines). Observational bias does not explain the flattening 
in Fig. 1 because the detection rate of radio afterglows is 
largely insensitive to redshift. Of the 60 GRBs with known 
redshift that have been observed in the radio, there are 42 
with detected afterglows. Their detection rate above and 
below z = 1 is identical. 

3. AFTERGLOW MODELING 

We now proceed to combine these radio data (§2) with 
the extensive X-ray and optical/NIR data for GRB 050904 
in order to constrain the physical parameters of the outflow 
and the circumburst medium. We interpret these multi- 
wavelength data within the framework of the relativistic 
blastwave model (see Meszaros 2002 for a review). The 
particular fitting approach that we take is described in 
more detail in Yost et al. (2003). 

12 The Very Large Array is operated by the National Radio Astronomy 
under cooperative agreement by Associated Universities, Inc. 



To limit the number of free parameters in the fit, we 
chose to model only the evolution of the forward shock 
after energy injection has ceased. The early X-ray light 
curve shows complex flaring activity, common in Swift 
bursts but lasting an unusually long time (At < 0.7 days), 
corresponding to a central engine lifetime of more than 
2 hrs in the source rest frame (Cusumano et al. 2006). 
On-going energy injection is also suspected based on the 
flares and the flattening of the optical/NIR lightcurves 
on similar timescalcs (Boer et al. 2006; Haislip et al. 
2006). If these early data are excluded the NIR data 
show a smooth power-law evolution with a clear break at 
tj = 2.6 ±1.0 days and decay indices before and after the 
break of ai = — 0.72±g;^ and a 2 = -2.4 ±0.4 (Tagliaferri 
et al. 2005). If interpreted as a jet break, the sharpness 
of the transition makes it unlikely that this burst occurred 
in a wind-blown environment (Kumar & Panaitescu 2000) 
and motivates our modeling choice of a constant density 
circumburst medium. 

Our best-fit forward shock model is shown in Fig. 2. We 
fit for the isotropic kinetic energy of the shock Ek.i SO ,52 , the 
opening angle of the jet 9j, the density of the circumburst 
medium n , the electron energy index p, and the fraction 
of the shock energy density in relativistic electrons e e and 
magnetic fields e^. 

Our model makes the unusual prediction that the syn- 
chrotron cooling frequency lies below both the X-ray and 
optical bands before the first J-band detection of the af- 
terglow at Ai=3.07 hrs (Haislip et al. 2006). This implies 
a steep spectral slope of (i = —p/2 = —1.07 (Sari, Piran & 
Halpcrn 1999). This is likely the case at At=3 days where 
we derive an NIR/X-ray slope of (3 ox = —1.1 ± 0.05. This 
value also agrees with Tagliaferri et al. (2005) who derive 
a best fit optical/NIR spectral index of (3 = —1.25 ± 0.25 
over the interval from 0.4 to 7 days postburst. A steep 
spectral slope of f3 x ~ —0.85 is seen even when the X-ray 
light curve is highly time variable between 0.05< At < 0.7 
days. Thus we find that the available data support a small 
value for the cooling frequency. 

Our best-fit model favors a circumburst density n o =680 
cm~ 3 . To test the robustness of this result we ran several 
model fits, fixing n over a range of values between 0.7 
cm~ 3 and 7xl0 4 cm~ 3 . As expected, the X-ray and NIR 
data are relatively insensitive to n , so the radio afterglow 
provides the tightest constraints. Large densities n > 10 4 
cm~ 3 are ruled out since the solutions are bounded by 
the condition that the magnetic and electron densities do 
not exceed equipartition. Low densities n < 10 cm~ 3 
are also not favored since they predict bright early radio 
emission in the first week well in excess of that which is 
observed (§2). Over the allowable density range, Ek,i SO ,52 
varies from 40 to 530, and therefore it is no better con- 
strained by the afterglow data than the isotropic radiated 
energy -E 7i i S0 .52 is by the prompt emission (Cusumano et 
al. 2006).' 

Another possibility is that radio afterglow flux density 
measurements are biased by the turbulent ionized inter- 
stellar medium of our Galaxy which induces non-Gaussian 
intensity fluctuations for compact radio sources. Indeed, 
it has been shown by Cordes & Lazio (1991) that repeated 

Observatory, a facility of the National Science Foundation operated 



short observations similar to those in Table 1 increase the 
odds of detecting a weak radio signal close to the detection 
threshold. We replaced the detections with upper limits 
and re-ran the fit. High density solutions are still preferred 
(n ~ 10 4 cm~ 3 ) because both the X-ray and optical after- 
glows of GRB 050904 were bright (Cusumano et al. 2006; 
Haislip et al. 2006). Thus, in order to suppress the ra- 
dio emission at 8.5 GHz (or 62 GHz in the rest frame), a 
large value of the synchrotron self-absorption is required 
and hence a large circumburst density. 

4. DISCUSSION 

With the detection of the radio afterglow from 
GRB 050904 we have a complete multi- wavelength dataset 
for this high redshift burst. By comparing the physical 
properties of the burst and its circumburst environment 
to those GRBs found at lower redshifts, we hope to gain 
some insight into the birth and death of the earliest gen- 
erations of stars in the Universe. 

There is good support for intcrprctating GRB 050904 
as an energetic event. From our afterglow modeling (§3) 
of the isotropic kinetic energy of the shock Ek y i SO ,h2 = 88 
(normalized in units of 10 52 erg). Estimates obtained for 
the isotropic radiated energy (Cusumano et al. 2006) with 
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the uncertainty in the location of the peak of the gamma- 
ray spectrum. These constraints can be refined owing to 
our fitting the jet opening angle 0j=O.14 rad (see also 
Tagliaferri et al. 2005). The geometrically corrected ener- 
gies are E k ^ 2 ^ E k , is0t52 x fl 2 /2=0.9, and 0.7 < £ 7)52 < 3. 
Both the radiated and kinetic energy of this event lie at 
or beyond the values derived for a large sample of lower 
redshift afterglows (Panaitescu & Kumar 2001, 2002; Yost 
et al. 2003; Friedman & Bloom 2005). An independent 
check on this result uses the X-ray light curve to com- 
pute the geometrically-corrected X-ray luminosity at some 
fiducial time (usually taken at At=10 hrs). We derive 
1^=2. 6x 10 45 erg s _1 , again larger than any previous GRB 
afterglow (Berger, Kulkarni & Frail 2003). 

Perhaps the most interesting aspect of these data is 
that they allow us to make an in situ measurement of 
the density structure for a massive star in the early Uni- 
verse. From our afterglow modeling we conclude that 



GRB 050904 exploded into a constant density medium 
with n o =680 cm -3 - two orders of magnitude above the 
nominal value for lower redshift GRBs (Friedman & Bloom 
2005; Sodcrbcrg et al. 2006). Line-of-sight measures 
(Cusumano et al. 2006; Haislip et al. 2006; Totani et al. 
2005) also indicate that a substantial column of gas exists 
towards GRB 050904. The fine-structure silicon lines Sill* 
are the most important since their column density ratios 
yield estimates of the electron density n e ~ io 2 - 3±0 - 7 cm~ 3 
(Kawai et al. 2006). The close agreement of this line-of- 
sight estimate with our in situ value argues for a common 
physical origin for the gas in the immediate vicinity of the 
GRB. 

There are testable predictions for star formation and 
collapse in the early Universe (Bromm & Loeb 2006 and 
reference therein). The first stars are expected to be very 
massive and could produce energetic GRB explosions like 
GRB 050904. Likewise, the high density we derived for 
GRB 050904 may be part of the general increase expected 
for the ambient density of the interstellar medium (Cia- 
rdi & Loeb 2000), but more likely just indicates that this 
GRB exploded in a dense molecular cloud. A lower density 
would be expected if the progenitor of GRB 050904 was a 
low metalicity Wolf-Rayet star (Vink & de Koter 2005). 
However, it remains a problem, at both low and high red- 
shifts, why the radial density signature from the collapsar 
wind is rarely detected (Chevalier, Li & Fransson 2004). 
A larger sample is needed before drawing too much from 
these predictions. Deep radio observations will continue 
to play a role since can be seen out to the highest red- 
shifts and are sensitive probes of the total calorimetry of 
the explosion and the density structure of the circumburst 
medium. Looking further ahead to the next generation 
of centimeter (EVLA) and millimeter radio instuments 
(ALMA), events like GRB 050904 will be bright enough to 
sample primordial magnetic fields (Gaensler, Beck & Fer- 
etti 2004), and probe the cold atomic and molecular gas 
from the time of the first light (Inouc, Omukai & Ciardi 
2005). 



DAF wishes to thank B. Clark for his generous alloca- 
tion of VLA time during the dynamic scheduling tests. 
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Table 1 
Radio Observations 



Date Obs 


At 


F v a 


(UT) 


(days) 


OJy) 


2005 Sep. 4.58 


0.50 


89±58 


2005 Sep. 5.48 


1.40 


41±25 


2005 Sep. 9.48 


5.40 


-3±25 


2005 Sept. 10.31 


6.23 


27±24 


2005 Sep. 24.20 


20.12 


89±37 


2005 Oct. 3.36 


29.28 


40±30 


2005 Oct. 7.30 


33.22 


-10±35 


2005 Oct. 8.26 


34.18 


64±23 


2005 Oct. 9.34 


35.26 


116±18 


2005 Oct. 11.43 


37.35 


67±17 


2005 Oct. 18.37 


44.29 


13±27 



a All errors are given as la (rms) 
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Fig. 1. — Radio flux density versus redshift for a complete sample of 60 bursts observed at radio wavelengths from 1997 to 2006. 
The radio measurements were made at 8.5 GHz and are taken from Frail et al. (2003) and the public radio afterglow database 
(http://www.aoc.nrao.edu/-~dfrail/grb_public.shtml). For clarity only the 42 GRBs with known redshifts and radio afterglow detections 
are plotted. The large (small) circles indicate measurements taken less (more) than 5 days before (after) the burst in its rest frame. The 
long dashed line is the flux density for a canonical afterglow model at At=5 days (in the observer frame). It assumes E^ i so 52 = 10 53 erg, 
9j = 0.1 rad, n = 10 cm -3 , p = 2.2, e e =0.1, and eb=1% (see text for more details). The short dashed line shows the expected decrease if 
the flux scales simply as the inverse square of the luminosity distance. 
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Fig. 2. — Broad-band afterglow fits for the radio, near infrared (NIR) and X-ray data of GRB 050904. The radio data are taken from §2. 
Solid circles indicate detections, upper limits are inverted triangles (plotted as flux + 2a). The open square is the peak flux obtained by 
averaging the first week of data. The NIR data are taken from Tagliaferri et al. (2005) and Haislip et al. (2006), with a small correction 
made for Galactic dust extinction (Schlegel, Finkbeiner & Davis 1998) before converting to flux density (Bessell, Castelli & Plez 1998). JHK 
and K' are shown for display purposes on the same plot after scaling by the best-fit spectral index. We have converted the Swift X-ray data 
( 0.2-10 keV) to flux density using the average photon index T = —1.84 and a frequency of 2.8xl0 17 Hz (see Cusumano et al. 2006). The 



cm , p=2.I4, e e =2.0%, and es=l-5% (see text for more details). The thin dashed line shows the effect of model fitting when the value of 
n is fixed at 100 times lower than the best-fit value. 



